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We report our observational results of 870 /im continuum emission and its 
linear polarization in the massive star formation site W51 e2/e8. Inferred from 
the linear polarization maps, the magnetic field in the plane of sky (Bj_) is traced 
with an angular resolution of 0'.'7 with the Submillimeter Array (SMA). Whereas 
previous BIMA observations with an angular resolution of 3" (0.1 pc) showed a 
uniform B field, our revealed B^ morphology is hourglass-like in the collapsing 
core near the Ultracompact H II region e2 and also possibly in e8. The de- 
crease in polarization near the continuum peak seen at lower angular resolution 
is apparently due to the more complex structures at smaller scales. In e2, the 
pinched direction of the hourglass-like B field morphology is parallel to the plane 
of the ionized accretion flow traced by H53a, suggesting that the massive stars 
are formed via processes similar to the low mass stars, i.e. accretion through 
a disk, except that the mass involved is much larger. Furthermore, our finding 
that the resolved collapsing cores in e2 and e8 lie within one subcritical 0.5 pc 
envelope supports the scenario of magnetic fragmentation via ambipolar diffu- 
sion. We therefore suggest that magnetic fields control the dynamical evolution 
of the envelope and cores in W51 e2 and e8. 

Subject headings: ISM: individual (W51 e2/e8) — ISM: magnetic fields — po- 
larization — stars: formation 



1. Introduction 



The magnetic (B) field has been suggested to play an important role in the star for- 
mation p rocess. While the B field flux density is eventually redistributed via ambipolar 
diffusion (IMestel fc Spitzer 19561 : iMouschovias 19781 ). the collapse itself is slowed sufficiently 
to explain the low star formation rate obsery ed in molecular clouds. Alternative support 
via turbulence (cf. iMac Low fc Klessen 2004 ) seems less important on parsec (pc) scales 
since the B fields in the plane of sky (B i ) are often obs erved to be orga nized and uniform 
across the cl oud, such as in M17 JDotson 1996h . OMC-1 JSchleuning 1998^ and DR21 MAIN 
(IKirby 20091 ). One key question is at which sizescale will the magnetic support be overcome 
by gravity. The morphology of the B field at that point may reveal the details of the con- 
traction process such as geometry and timescale. The Submillimeter Array (SMA) can be 
used to address this question by resolving the B^ structures via dust polarization studies 
with high angular resolutions at typically a few arcseconds. 



The B field is traced by the dust continuum emission. The dust grains are most likely 
not spherical in shape, but somewhat elongated. They are thought to be aligned with their 
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minor axes parallel to the B field in most of the cases (ILazarian 20071 ). Among different 
alignment mechanisms, radiat ion torques seem to be a pr o mising mechanism to al ign the 
dust grains with the B field (IDraine fc Weingartner 19961 ; ILazarian fc Hoang 20071 ). Due 
to the differences in emissivity perpendicular and parallel to the direction of alignment, 
the observed thermal dust emission will be linearly polarized. The direction of the linear 
polarization is therefore perpendicular to the B field. With the SMA, we are able to detect the 
polarized component of the thermal dust emission at sub-millimeter (sub-mm) wavelengths 
in order to trace the B field within the dense cores, where stars are formed. Compared 
to the polariz ation studies via abso rption and scattering of stellar light in the optical or 
near infrared (IGoodman et al. 19951 ). sub-mm polarization, being derived directly from dust 
emission, does not suffer from a limited range in grain size and the possible contamination 
from the more diffuse emission and absorption along the line of sight. 

In this paper, we present SMA observational results with an angular resolution of 0'.'7 
(0.02 pc) of the massive star forming site W51 e2 and e8 in W51 MAIN. The dust continuum 
at a wavelength of 870 /im, and the B^ field inferred from its linearly polarized component 
are p resented. The W51 MAIN is on the eastern edge of W 51. It is at a distance of 7. Oil. 5 
kpc (IGenzel et al. 19811 ) or 6.1±1.3 kpc (jlmai et al. 20021 ). Here, we adopt a distance of 7 
kpc. There is a group of Ultracompact H II (UC HII) regions in W51 MAIN, and many H ? Q , 
OH and NHs mase r spots have been identified (IGenzel et al. 198ll ; iGaume &: Mutel 19871 ; 
Pratap et al. 19911 1 to be associated with the e2 and e8 regions. The terminology of the struc- 
tures discussed in this paper is shown in the schematics in Figure 1. The radio continuum 
sources e2, e4, e8, el and e3 are UCHII regions (IGaume &; Johnston 19931 ; IZhang Ho 19971 ). 
and their locations are labelled in Figure 2. Hereafter, e2, e4, e8, el and e3 (when 
in italic) refer to their corresponding UCHII regions. The infall signatures toward the 
e2 and e8 regions (i.e., e2 and e8 c ollapsing cor es) have been det ected clearly in NH 3 
JHo fc Young 19961 : Izhang fc Ho 1997h and in CS Jzhang et al. 1998T ). indicating that they 
are in an e arly evolutionary stage. Furthermore, the total luminosity of the W51 MAIN is 
2xl0 6 L (jJaffe et al. 19871 ). indicating that it is a massive star forming site. 



The polarized dust emissions associated with the envelope of the W51 e2 and e8 regions 
have been previously observed at 1.3 mm and 850 /im. The B^ field structure varies with 
different size scales. Chrysostomou et al. (2002) has shown that the morphology of the field 
on the very large scale observed with SCUBA with an angular resolution of ~ 10" (0.5 pc) 
appears more complex, possibly because of projection effects from several clouds along the 
line of sight. With an angular resolution of ~3" (0.1 pc) with the BIMA, Lai et al. (2001) 
found that the position angles (P.A.s) of the polarization vectors vary smoothly across the e2 
and e8 cores (Figure 2(a)), suggesting that the B field dominates over the turbulent motions 
in the envelope. At which scales will the B field lose its dominance over turbulence and 
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gravity? 



Observation 



The observations were carried out on 2008 July 13 using the SMA ( Ho et al. 2004 1"*] in 



the extended configuration, with seven of the eight antennas available. The projected lengths 
of baselines ranged from 30 to 262 kA. The largest size scale which could be sampled in this 
observation was ~8" (0.3 pc). The local oscillator frequency was tuned to 341.482 GHz. 
With the 2 GHz bandwidth in each sideband, we were able to cover the frequency ranging 
from 345.5 to 347.5 GHz and from 335.5 to 337.5 GHz in the upper and lower sidebands, 
respectively. The phase center is near e2 at Right Ascension (J2000) = 19 h 23 m 43 s .95, 
Declination (J2000) = 14 o 30'34'.'00. e8 is ~7" south of the phase center. The primary beam 
(field of view) of the SMA at 345 GHz is ~30". 

Linear polarization (LP) observations using interferometer arrays are best obtained 
using receivers which detect both orthogonal circular polarizations (CP) simultaneously. 
However, the SMA receivers are intrinsically linearly polarized and only one polarization is 
available currently. Thus, quarter-wave plates were installed in order to convert the LP to 
CP. Detailed information of the design of the quarter-wave plates and how the quarter-wave 
plates were controlled is described in Marrone et al. (2006) and Marrone & Rao (2008). We 
assume that the smearing due to the change of the P.A.s on the time scale of 5 minutes in 
one cycle of polarization measurement is negligible. 

The conversion of th e LP to CP is not perfect. This instrumental polarization (also called 



the leakage terms) (see ISault et al. 19961 ) and the bandpass were calibrated by observing 
3c454.3 for 2 hours while it was transiting in order to get the best coverage of parallactic 
angles. The instrumental polarization is ~1% for the upper sideband and ~3% for the lower 
sideband before calibration, and ~0.6% after calibration in both sidebands. The complex 
gains were calibrated every 12 minutes by observing 1751+096 until it set, followed by 
1925+211 for the last 3.5 hours. The absolute flux scale was calibrated using Titan. 

The data were calibrated and analyzed using the MIRIAD package. After the standard 
gain calibration, self-calibration was also performed by selecting the visibilities with uv 
distances longer than 40 kA. In order to Fourier transform the measured visibilities to the 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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image, the task INVERT in MIRIAD was used with natural weighting. The Stokes Q and U 
maps are crucial for the derivation of the polarization. We use the dirty maps of Q and U to 
derive the polarization in order to avoid a possible bias introduced from the CLEAN process. 
We applied CLEAN to the Stokes I (total intensity) map in order to reduce the sidelobes. 
The presented SMA images have all been corrected for the primary beam attenuation. The 
synthesized beam of the presented maps is &!7x&!6 with a P. A. of —58°. The presented 
polarization vectors are gridded to a 0'.'3 spacing - which is about half of the synthesized 
beam FWHM - in order to show the curvature of the B field morphology. Therefore, adjacent 
polarization vectors are not formally independent within one synthesized beam. However, 
as usual, relative information can be extracted at under the synthesized beam resolution. 

The Stokes /, Q and U images of the continuum are constructed with natural weighting 
in order to get a better S/N ratio for the polarization. The noise levels of the /, Q and U 
images are ~ 60, 4 and 4 mJy Beam -1 , respectively. The strength (J p ) and percentage (P(%)) 
of the linearly polarized emission are calculated from: 1^ = Q 2 ' + U 2 — Cq^u and P(%) = 
respectively. The term oqjj is the noise level of the Stokes Q and U images, and it is the 
bias correction due to the positive measure of I p (Leahy 1989; Wardle & Kronberg 1974). 
The <jj p is thus 4 mJy beam -1 . To derive the polarization, the MIRIAD task IMPOL was 
used. The SMA polarization vectors presented are above 3<7/ p in red segments and between 
2 to 3 07 p in black segments. 



3. Results 



The 870 jum continuum emission and its polarized components were detected (Figure 2; 
Table 1 & 2). The results are presented in this section. 



3.1. Continuum Emission 



In e2, a compact 870 //m continuum emission structure with a radius of ~ 1" (0.03 pc) is 
centered at ~0'. / 7 east of e2. Extending to the north-west of this compac t emission, a fainter 
structure with an over all length of ~ 2" (0 .07 pc) is detected. The H 2 (IGenzel et al. 19811 ) 
and (J,K) = (9,6) NH3 (IPratap et al. 199ll ) masers are located in this north-west extension, 
~ 2" away from the continuum peak. Associated with the continu um peak, there are OH 
masers detected within 0'.'5 to the east and ~1" to the south of e2 f lGaume fc Mutel 19871 ? 
), suggesting that it is an active star forming site. 



In e8, the 870 /im continuum peak is centered at 0'.'3 west of e8. e4, el and e3 are at the 
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periphery of the 870 /tin continuum emission. There is an extension toward the south-west 
with an overall length of ~3". Associated with e8, an NH 3 maser spot was detected 0'/8 south 
of the 870 /im continuum pe ak by Pratap et al. (1991). The OH (IGaume fc Mutel 19871 ) 
and H 2 (IGenzel et al. 19811 ) masers are also associated with e8 and the 870 /mi continuum 
peak, suggesting again that this is an active star formation site. 

When fitted with a Gaussian, the deconvolved size of the 870 /im emission in e2 is 
ff!9x0'!8, slightly larger than the synthesized beam, and therefore, the e2 core has been 
resolved. For e8, the deconvolved size is (f!9x0'!3 with a P.A. of 12°. Therefore, e8 has been 
resolved along the major axis of the dust ridge but not along the minor a xis. In both e2 



and e8, the 870 /a n continuum emissions are associated with the NH3 cores (IHo et al. 1983 
Zhang & Ho 19971 ). suggesting that they are also tracing the dense regions. 



The measured 870 /im flux densities within the upper and lower boxes in Figure 2(b), 
associated with e2 and e8, are 9.3 and 4.0 Jy, respect ively. The flux densities o f the free-free 
continuum Fff at 1 .3 cm in e2 and e8 are 300 mJy (IGaume fc Johnston 19931 ) and 17 mJy 
(IZhang fc Ho 19971 ). respectively. In order to estimate the Fff contribution at 870 /im, we 
extrapolate from 1.3 cm, assuming Fff oc v~ 01 . Although this assumption of optically thin 
emission is crude, it has been shown that the resultant Fff roughly agr ees (within a facto r 
of 3) with the estimate from the radio recombination line at 2 mm fjZhang et al. 19981 ). 
suggesting that the assumed Fff oc v~ 0A is reasonable. The extrapolated Fff at 870 /tm is 
~230 and 13 mJy for e2 and e8, respectively. As compared to the 870 /im flux densities, 
Fff contributes ~2% for the e2 region and 0.3% for the e8 region. Therefore, the 870 /tm 
continuum is dominated by dust emission. Hereafter, the structures traced by the 870 /tin 
emission in the e2 and e8 regions are named as e2 dust ridge and e8 dust ridge, respectively. 



Assuming a dust temperature of 100 K (IZhang et al. 19981 ) . a dust grain emissivity 
Q(X) oc X' 13 with (3 = 1, and the normal gas to dust ratio of 100, we estimate gas masses 
M gas of 245 and 106 M Q for the e2 and e8 dust ridges, respectively (cf. Tang et al. 2009). 
Note that the M gas given here is highly affected by the assumed (3. If the assumed (3 is 2, 
the estimated M gas will be 14 times larger. Assuming the extents along the line of sight 
are equal to the diameters of the emission area in the e2 and e8 dust ridges, the average 
gas number densities nn 2 are 3.4xl0 6 and 2.2xl0 6 cm -3 , respectively. By using the same 
equation and the same ass umed values of (3 a nd dust temperature, the M gas estimated from 



the 2 mm dust continuum (IZhang et al. 19981 ) for the e2 and e8 dust ridges are 1100 and 590 
M Q , respectively. The difference in the estimated M gas at 2 mm and 870 /im is most likely 
due to the missing flux from the extended component, which is not recovered with our SMA 
observations. In comparison, with t he same assump tions, the M gas of the envelope is 1834 



M Q as traced at 1.3 mm by BIMA (ILai et al. 20011 ). The M gas associated with the e2 and 
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e8 dust ridges recovered with our SMA observations is ~ 19% of the M gas in the envelope. 

The main conclusion from the dust continuum data is that the associated mass is large. 
The morphology of the dust continuum is elongated. The positional offsets between the 
various embedded sources are significant, such as between the positions of the 870 //m peaks 
and the UCHII regions. These results are consistent with the formation of a cluster of stars. 



3.2. Dust Polarization 

The polarization in the e2 and e8 dust ridges is detected and resolved (Figure 2 (c) and 
(d)). Throughout the paper, P. A. is defined from the north to the east. In the e2 dust ridge, 
the bulk of the polarization vectors form a ring around the 870 /im peak with a radius of 
~1" and with the geometric center near the continuum peak instead of e2. In the north-west 
extension of the dust ridge, the polarization appears to be perpendicular to the major axis 
of the extension. 

The e8 dust ridge is ~7" away from the phase center. Even though the antenna response 
is 15% less efficient than at the phase center, the polarization revealed is clearly also not 
as uniform as previously seen with BIMA. The polarization vectors again form a ring like 
structure around the continuum peak. The polarization is weaker in e8 with more vectors 
between 2 to 3 oj p . 

In comparison, the polarization in the envelope of the e2 and e8 regions, as revealed 
with an angular resolution of 3" (0.1 pc) with BIMA, shows a relatively uni form distributio n 



in P.A. and therefore, a fairly uniform B^ field at 1.3 mm (Figure 2(a); lLai et al. 20011 ). 
In their results, the polarization in the e2 region is weak and resolved into e2 main and e2 
pol NW, named in the same paper, according to the P.A. of the polarization vectors. The 
component e2 pol NW is at 3" to the north-west of e2. There is a gap where no polarized 
emission is detected between e2 and e2 pol NW. In the e8 region, the polarization in the 
BIMA results is nearly uniform with a decrease in polarization percentage near the peak 
position. 

In order to test if the differences in polarization properties from SMA and BIMA are due 
to their different angular resolutions, we smoothed our SMA results to the BIMA resolution, 
as shown in Figure 3. Wherever the polarized emission was both detected at 1.3 mm and 
870 /im, the resultant P.A.s of the polarization differed by ~ 30° on average. This significant 
difference can be due to the different sampling of the visibilities, which are in the range of 6 
to 170 kA (A=1.3 mm) for the BIMA and in the range of 30 to 262 kA (A=870 /im) for the 
SMA. Specifically, the SMA filtered out the more extended and uniform component which is 
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larger than 8". At the same angular resolution, the derived global field directions in e2 
and in e8 are therefore consistent in the regions where both the SMA and BIMA have polar- 
ization detections. Most importantly, the smoothed SMA polarization map shows that the 
polarization percentage has decreased significantly, especially near the continuum peak posi- 
tions, where the field geometry is more complex at the resolution of 0'/7. This demonstrates 
that the low polarization percentage at the emission peaks is due to the limited angular 
resolution when a more complex underlying B field morphology has not been resolved. This 
effect can also be due to the decrease of the alignment efficiency of the dust grains in denser 
regions (Lazarian & Hoang 2007) or due to geometrical effects, such as the differences in the 
viewing angles (Gongalves et al. 2005). However, in this case, the complex B field structure 
is the dominant effect. 

The polarization percentage P(%) decreases with increasing continuum intensity I in 
both e2 and e8 even for the higher resolution SMA results (Figure 4). Since the BIMA results 
come from a resolution effect, the same might be true for the SMA results at the emission 
peaks. Away from the emission peaks, the general increase in P(%) is somewhat misleading. 
Figure 2(a) shows that this effect is not symmetrical on either side of the elongated envelope, 
i.e. the P(%) differs with positions on the same contour level of I. This is reflected by the 
large dispersion in P(%) at any value of I / I max . Several effects, including B field geometry 
related to the line of sight, need to be disentangled. That the P{%) ranges mainly between 
1% to 10% (Figure 4), seems to agree with the model of grain growt h in the dense region s 



where grain alignments are via radiative torques (see Figure 11 in iPelkonen et al. 20091 ). 
However, based on our results, the effects of angular resolution and geometry must first be 
taken into account. 



4. Discussion 

4.1. Hourglass B field Morphology inside the e2 dust ridge? 

The inferred B field in the e2 dust ridge exhibits a complex but organized morphology 
(Figure 5). We have tested the hypothesis of the measured B field being radial, and have 
shown quantitatively the preference of a non-radial field at a high significance level (see 
Appendix). There are positions where no polarized emission (depolarization) is detected, 
extending along a P. A. ~ 60° across the 870 fim peak (color scale in Figure 2(c)). The 
existence of non-radial field lines together with the depolarized zones are in favor of an 
hourglass field morphology. Along the extension of the dust ridge toward the north-west, 
the B field lines are approximately parallel to the major axis, which is consistent with the 
BIMA measurement at 1.3 mm. These lines are radial-like, but the complex structure in the 
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north-west could belong to another embedded source as possibly indicated by the masers. 



Associated with e2, organi zed motions in the ioni zed gas have been revealed with the 
H53a radio recombination line (IKeto fc Klaassen 20081 ). with the maximum velocity gradi- 
ent along P. A. ~ 60°. These authors interpret this gradient as a supporting evidence for 
an accretion flow along a dense flattened structure, where the detected motion tracks the 
ionized particles on the surface of the dense midpla ne. Both the infal l and rotation near 
e2 have also been detected in several molecular lines (IHo fc Young 19961 ; IZhang &: Ho 19971 ; 
Zhang et al. 19981 ). As discussed in Keto & Klaassen (2008), this H53a accretion flow in 
the direction of PA. ~ 60° might drive the molecular outflow at P.A. ~ — 20° as traced 
by the CO 2-1 line. The argument is based on the hypothesis that if the massive star 
formation process is similar to the low mass case, the bipolar outflow should be along the 
rotation axis. The linearly dis tributed H 2 Q and OH masers in the W51 e2 region could 
trace an outflow (Figure 19(c), iDe Buizer et al. 20051 ). as identified with the CO 2-1 line. 
Although the determined direction may be highl y uncertain, the ro tation in NH3 (3,3) is 
more clear ly revealed along P A=135° (Figure 7 in IZhang fc Ho 19971 ) and in CH3CN along 
PA=110° (IZhang et al. 19981 ). which seems inconsistent with the gradient detected with the 
H53a ionized flow. This might indicate that the revealed kinematics based on different lines 
may be from multiple embedded sources. Higher spatial resolution kinematic studies with 
hot core molecular lines will be helpful for deciphering the underlying structures. 

The B field appears to be hourglass-like near e2, with the field lines pinched along 
the plane of the proposed H53a accretion disk. If the B field lines are frozen into the 
ionized material, the field lines will be tangled along with the rotation and infall motions. 
The revealed depolarization might then result from the more complex underlying B field. 
We note that the field lines seem to go to the core with an essentially radial pattern, and 
therefore, leading t o a sharp pinched angle in the hourglass. In contrast, the low mass case 
(IGirart et al. 20061 ) shows a wider and smoother pi nched angle. We sp eculate that a larger 
infall momentum and a larger differential rotation (IZhang et al. 19981 ) might drag the field 
lines along and result in a narrower pinched angle in the projected plane. Projection of a 
nearly pole-on hourglass-like morphology possibly also leads to similar signatures. In any 
case, the scenario of material accreting through a disk as proposed by Keto and Klaassen 
(2008) is supported by our inferred B field morphology. 



4.2. Hourglass B field Morphology in the e8 dust ridge? 



Along the e8 dust ridge, the B field also shows a systematic deviation from the larger 
scale (0.5 pc) B field revealed by BIMA. This can again be explained by the field lines being 
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dragged along with the accretion toward e8. In this case, the revealed B field appears to 
be part of an hourglass structure on a larger scale of 4" (~0.08 pc) (Figure 5(b)), with its 
pinched direction parallel to the dust ridge. Centered on the e8 continuum peak, a compact 
hourglass structure would be more convincing except for the field lines to the north. There 
are H2O masers north of the e8 continuum peak, and another embedded source may be 
indicated. This could explain the incomplete hourglass structure here. 

As in the case of e2, a zone of depolarization seems to be present at the continuum 
peak, along the north-south direction. This is consistent with the pinch direction of the 
hourglass-like morphology being along the elongated e8 dust ridge. Rotation associated 
with the e8 collapsi ng core was detected in the direction of P.A. ~ 156° (—24°) with CH3CN 



(IZhang et al. 19981 ). In this scenario, the pinch direction of the hourglass-like B field is 
parallel to the plane of rotation. The rotation axis of the e8 collapsing core is then almost 
parallel to the B field threading the 870 (im dust ridge . Note that the rotation direction as 



traced in CH3CN is still uncertain (IZhang et al. 19981 ). An accurate determination of the 



plane of rotation associated with the e8 collapsing core is needed to test if the larger scale 
B field controls the direction of accretion. 

Although the plane of accretion (or the pinched angle of the ho urglass) cannot be 



determined with certainty, the collapse signature was detected toward e8 ( IHo fc Young 1996 



Zhang fe Ho 19971 ; IZhang et al. 19981 ). consistent with the possible hourglass-like B field 
morphology. Furthermore, this collapsing core is inside the 0.08 pc scale dust ridge, as 
revealed with the C/!7 angular resolution B field morphology. Based on this morphology 
and the presence of e4, e8, el, and e3, we suggest that the star formation process involves 
different stages of fragmentation, proceeding at different evolutionary timescales. 



4.3. Estimate of the Strength of the B field 

The B field strength can be estimated by comparing the gravitational force /g with 
the B field tension f B following Dotson (1996) and Schleuning (1998). The value of /q at a 
distance Rq away from the center is given by 

_ GM R p _ 26 M R n U2 R G 2 dyne 

/G R 2 G lOOM lO 5 cm- 3l O.lpc J cm 3 ' [) 

where Mr refers to the gas mass enclosed within a radius Rq, p is the mass density at Rq, 
and 77,h 2 is the gas volume number density. The /b can be given by 

1 -> - B 2 9K/ B , 9 . R B . , dyne . . 

/ B = j-B ■ yB ~ — — = 5 x Kr 26 (— ) 2 (— ^-J' 1 (2) 



Air ' AttR-q mG 0.5 pc cm 
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where Rb is the radius of a magnetic flux tube, and B is the B field strength. Since the 
e2 and e8 cores are known to be in a collapse stage, we conclude that / G > / B . An upper 
limit of B can then be derived. For the e2 collapsing core, Mr is estimated to be 22OM 
based on the 870 /im flux density within a radius of 1" of t he continuum peak . This is 
consistent with the proposed self gravitating mass of >160 M (|Ho fc Young 19961 ) based on 
the kinematics of the NH3 lines. Rq is 1" (0.034 pc), and the mean nn 2 within Rq is 2.7xl0 7 
cm -3 . Assuming Rq = R^ — 0.034 pc, the B field strength in the e2 core is therefore < 19 
niG. In the e8 collapsing core, Mr is 94 M Q , and % 2 is 1.2 xlO 7 cm -3 within a radius of 
1" centered on the peak position. The B field strength in the e8 core is therefore < 8 mG. 
Both of the upper limits of B are consistent with the lower limit of the larger scale B^ field 
strength of 1 mG (jLai et al. 2001 ) estimated from the method proposed by Chandrasekhar 
& Fermi (1953). 



4.4. Characteristic Length Scales 

To analyze the interactions between B field, gravitational force and thermal force in star 
forming sites, we further calculate the following three length scales following Mouschovias 
(1991): First, the interplay between ambipolar diffusion and Alfven waves is characterized 
by the Alfven length scale Aa- Second, the interplay between gravitational and thermal 
pressure forces is characterized by the thermal Jeans length scale At^d following Bonnor 
(1956) and Ebert (1955; 1957). Third, the interplay between magnetic and gravitational 
forces is characterized by the critical magnetic length scale Am,ct- They can be calculated 
using the following equations: 



and 

A M ,cr = 36— p7( in6 nH2 -s )" 1 m P c - ( 5 ) 
mG 10° cm 6 

Here, the parameter k (Eq. (6f) in Mouschovias 1991), related to the mean collision time 
between an ionized and a neutral particle, is assumed to be 0.5 when we derive Eq. (3), 
which is within the most likely range given in the reference in their paper. The f actor K is 



related to the cosmic ionization rate. We assume K = 3xl0~ 3 (IMouschovias 19911 ) following 
the ionization rate calculated by Nakano (1979). With the assumed k and K, the estimated 
fractional ionization rate is 3xl0~ 9 for a number density of 10 7 cm -3 , which seems to be 
reasonable. T is the gas temperature, and nn 2 is the gas volume number density. B is the B 
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field strength. Note that T is assumed to be 100 K in both the e2 and e8 dust ridges based 
on the analysis of the hot core lines by Zhang, Ho & Ohashi (1997). Since these natural 
length scales depend on n Ha and B, they are calculated separately based on the detected 
continuum emission with the same assumption as in §3.1. In the 1.3 mm envelope, the 
is the mean number density within a best-fit Gaussian centered on the peak, and B is the 
lower limit of 1 mG. In the 870 fim dust ridges, nn 2 is calculated within a radius of 1", and 
B is the upper limit calculated in §4.3. The calculated natural length scales in the e2 and 
e8 cores are listed in Table 3. 

The physical meaning of these length scales are explained clearly in Mouschovias (1991) 
and references therein. Aa gives the lower limit of the scale at which the B field can sustain 
the structure. At the scale R< Aa, the ambipolar diffusion between neutral and ionized 
particles is more efficient and the Alfvenic motion is less important. At, ct gives the scale 
where the gravitational force is equal to the thermal pressure. If an object has a size scale 
R- > AT, cr , gravity can overwhelm the thermal pressure and collapse will start. Am.ct gives 
the upper limit of the scale where the cloud can be magnetically supported along the B field 
direction. In a region R> Am,ct> there will be enough mass and therefore the material can 
collapse if there are no other supporting forces. 



4.4.I. Correlation with the SMA e2 dust ridge 

In the e2 dust ridge, the revealed B^ morphology is clearly pinched with a radius 
of ~ 0'.'8 near e2, comparable to the radiu s of the proposed rotating flattened structure 
of 1 " Jzhang fc Ho 1997 ; Zhang et al. 1998 ) and the proposed ionized accretion disk of ~ 



C'5 (IKeto fc Klaassen 20081 ). The derived AT, cr is ~ 0'/2, suggesting that at ~ 1", gravity 
will easily overcome the thermal pressure support if there are no other supporting forces. 
The scale where ambipolar diffusion starts to take place (Aa) is ~ 0'.'2 for the e2 dust 
ridge, consistent with the observed pinched B^ field lines. Note that the revealed width of 
the depolarization zone near e2 is narrow (< 0'.'5), which is smaller than our synthesized 
beam. Higher angular resolution measurements with at least 0'.'3 resolution are needed 
to discriminate whether the depolarization is due to ambipolar diffusion, inefficient grain 
alignment or other mechanisms, such as geometrical effects. The calculated scale where the 
B field can sustain the structure against gravitational collapsing (Am, ct ) is 0'.'8 along the field 
line. However, it is difficult to compare Am.ct with the scale associated with the dust ridge, 
because the large scale (0.5 pc) B^ field is twisted by ~45° (Figure 2(a)) at 3" resolution. 
Observations with visibilities at both shorter and longer uv ranges are needed in order to 
link the B^ in the core with the field in the envelope at the same wavelength. Note that the 
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weak constraints on Aa and Am,ct result from the large range of possible B field values. 



4-4- 2- Correlation with the SMA e8 dust ridge 

The dust continuum emission appears to be ridge-like, and the minor axis of the e8 dust 
ridge is approximately parallel to the 0.5 pc field direction. The deconvolved length of 
the e8 dust ridge along the minor axis is barely resolved, and we adopt an upper limit to the 
radius of 0'.'3 along the minor axis. This is consistent with the estimated Am, ct < CK'7, and 
the estimated AT, cr of 0'.'3. This suggests that thermal pressure is significant as compared 
to gravity and field tension at this scale. Along the major axis, the deconvolved size is 0'.'9, 
which is larger than At, ct - Furthermore, ambipolar diffusion is expected to dominate at the 
scale Aa < 0'.'2 in e8. Hence we expect collapse and fragmentation to occur along the ridge. 
This is consistent with the revealed hourglass-like B field morphology associated with the e8 
collapsing core at &!7 resolution and the smooth B field morphology in the envelope at 3" 
resoluti on. These results in the e8 co re seem to be consistent with the ambipolar diffusion 



model (IMouschovias &: Morton 199ll ) and suggest that the formation of the dust ridge is 



influenced by the B field in the envelope. 



4-4-3 ■ Connection between the BIMA 0.5 pc scale envelope and the SMA dust ridges 

The main difference between the 0.5 pc envelope and the dust ridges is in the large 
contrast in the derived value for % 2 . Nevertheless, the agreement in the estimates of Aa 
and Am,ct for the envelope and the ridge is very good, because of the lower estimates of B in 
the envelope. This implies that the B field support is adequate till the scale of C'2— 1" . The 
derived value of At, ct is larger for the envelope than for the ridge, but still smaller than the 
measured size of the envelope. This is also consistent with the proposed scenario that the B 
field provides the support against gravity at the 0.5 pc scale. 



4.5. Role of Bi from Envelope (0.5 pc) to Collapsing Cores (0.02 pc) 

The structure of the B field varies with different size scales (Figure 6). 

At the 0'.'7 (0.02 pc) scale, the B field in the W51 e2 and e8 cores is not uniform. In 
both the e2 and e8 cores, the hourglass-like morphology suggests that at the 0.02 pc scale, 
gravity dominates over the B field. Our estimated collapsing M gas is on the order of 100 M 
(see §4.3) in both the e2 and e8 cores, which is roughly consistent with M gas estimated in 
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the local collapsing scenario (IHo &: Young 19961 ; ISollins et al. 20041 ). The cores are therefore 
in the supercritical phase. 

At the 0.5 pc scale, the B field in the envelope is uniform throughout the e2 and e8 dust 
ridges, except to the north-west of e2. Due to the small dispersion of the measured P.A.s 
of the polarization in the envelope, the B field is suggested to dominate over the turbulence 
with a strength of >1 mG by Lai et al. (2001). We further suggest that the envelope is 
subcritical. The reasons are: 1. B field dominates over turbulence in the envelope. 2. B field 
is apparently sufficient to support the envelope against gravity (§4.4.3). 3. The collapse is 
apparent only locally in e2 and e8 with M gas in the order of 100 M Q , so that the envelope is 
stable. 



Based on the MHD simulations (IKlessen et al. 20001 ). turbulent motions (cf. Mac Low & 
Klessen 2004) can produce an elongated envelope and can sustain the envelope from collapse. 
However, such elong ated structures would not exhibit a preferred alignment with the B field 
( jHeitsch et al. 20011 ). as it is seen here. The measured B field structures from 0.02 pc to 0.5 
pc do not show an obvious necessity for turbulent support. Instead, the cloud morphology, 
sizescale, and B field geometry are consistent with magnetic fragmentation via ambipolar 
diffusion jMouschovias fc Morton 199 it iLizano fc Shu 19891 : Ishu et al. 20041 ) . 



4.6. Comparison with other star formation sites 

Because the collapse signatures in the W51 e2 and e8 regions are clearly detected, and 
the UCHII regions are still relatively compact and weak, they are in an earlier evolutionary 
stage as compared with the other massive star formation sites, such as G5. 89— 0.39 (Tang et 
al. 2009). In G5. 89-0. 39, the B field is suggested to be overwhelmed by the turbulent motions 
from the UCHII expansion and the molecular outflows. A more complex B^ morphology is 
detected (Tang et al. 2009) with a spatial resolution of 0.02 pc. In contrast to G5. 89— 0.39, 
the B± morphologies are hourglass-like in both of the W51 e2 and e8 regions at the 0.02 pc 
spatial resolution, but much smoother at the 0.5 pc scale. This comparison of G5. 89-0. 39 
and W51 e2/e8 with the same spatial resolution indicates that the role of the B field varies 
with the evolutionary stages of the central sources. 

In the low mass star formation region NGC 1333 IRAS 4A, the hourglass shape is 
observed on a scale of 2400 AU (0.01 pc) and M gas is ~1.2M Q (Girart, Rao, & Marrone 
2006). In comparison, the hourglass B^ structure detected in W51 e2 is on the scale of ~ 
0.03 pc (1"), and the mass involved is ~ 200 M Q . The consistency of the directions of the 
ionized flow and the pinched field further suggests that the stars are formed with similar 
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mechanisms, i.e. material is accreted through a flattened structure. The difference is that 
the scale and the mass involved are much larger in the massive star forming regions. At 
the time of publication of this paper, additional observational evidence of an hourglass B 
field morphology in the massive star forming core G31.41+0.31 is presented in Girart et al. 
(2009). This further supports the proposed similar formation mechanism as in the low mass 
case. 

However, the massive star forming sites are much further away than the low mass regions 
in general. For example, W51 e2 is at 7 kpc, ~ 23 times further away than NGC 1333 IRAS 
4A. Observations with higher spatial resolution are thus needed for the massive star forming 
regions. The closest massive star forming site is Orion BN/KL. Source I in Orion BN/KL 
is suggested to be in an early stage of ma ssive star formation due to the weak and compact 



free- free emissi on (IPlambeck et al. 19951 ). Source I is also suggested to harbor an ionized 



accretion disk ( Reid et al. 2007h . The B^ at 0.5 pc starts to exhibit a larger scale hourglass- 



like morphology (jSchleuning 19981 ). The existence of the ionized disk, the uniform large scale 
B_|_ geometry and the compact free-free source all suggest that it is an analog of W51 e2/e8, 
but at a much closer distance. The comparison between w51 e2/e8 and Source I in Orion 
BN/KL may provide a clue of the B^ morphology at even higher physical resolutions. 

In W51 e2/e8, we found two supercritical cores at 0.03 pc within a subcritical envelope 
at 0.5 pc (§4.5). At a larger scale, DR21 MAIN is su ggested to h ave started undergoing a 



gravitational collapse in the central part of the cloud (IKirby 20091 ). Suggested by the same 
author, in the outer part (~ 1 pc away from the center), the cloud is still magnetically 
supported. That the collapsing cores formed inside a magnetically supported cloud as in 
W51 e2/e8, is therefore not a special case. As indicated by Vaillancourt (2009), observations 
of B field at the larger scales are needed to test the magnetically controlled star formation 
process. Shorter spacing visibilities of the SMA are needed to directly compare the field 
morphology of the cores with the envelope at the same wavelength of 870 fim. In this paper, 
we have demonstrated that smaller scale B^ maps provide crucial information about the 
B± field by resolving the star forming cores and by linking the field morphology with the 
kinematics of the molecular cloud. 



5. Conclusion and Summary 

To study the role of the B field in the star forming cores, we have observed and analyzed 
the B^ morphologies, inferred from the linearly polarized dust continuum emission, in the 
massive star forming site W51 e2/e8 by using the Submillimeter Array (SMA). We further 
compare the B^ morphologies in the dust ridges with the one in the envelope. Three different 
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natural length scales, namely the Alfven length scale Aa, thermal Jeans lengths scale Ax,cr, 
and the magnetic length scale Am,cd are calculated and compared to the dust ridges and 
envelope. Here are the summary and conclusion. 

1. The 870 /zm continuum and its polarization in W51 e2 and e8 are resolved with an 
angular resolution of ft' 7 as observed with the SMA. The polarization in both e2 and 
e8 exhibit complex structures. In comparison, the polarization at 1.3 mm observed at 
3" resolution with the BIMA revealed a uniform morphology a cross e2 and e8 , with 



almost no polarization detected near the peak position in e2 (lLai et al. 20011 ). We 
conclude that low or no polarization near the emission peaks heretofore seen in the 
star formation regions is likely due to a more complex underlying morphology 
(§3-2). 

2. In the e2 dust ridge, the inferred B^ morphology is hourglass-like near the collapsing 
core, with its pin ched direction parallel t o the direction of the ionized accretion flow as 



traced by H53a (IKeto fc Klaassen 20081 ). The B^ here shows a similar morphology as 
in the low mass star formation case NGC 1333 IRAS 4A. However, the mass included 
in this core is ~200 times larger. This result shows that the B field in the e2 collapsing 
core plays a similar role as in the low mass star formation regions (§4.1). Higher 
angular resolution observations are required to test if the hourglass-like B^ morphology 
is preserved at smaller scale. 

3. The e8 dust ridge is perpendicular to the 0.5 pc scale B^ direction as revealed with the 
BIMA, which suggests that the B field at 0.5 pc controls the forming process of the 
dust ridge. The B^ along the dust ridge exhibits a systematic deviation from the B^ 
at 0.5 pc scale. Associated with the e8 collapsing core, the hourglass B^ morphology 
is more clearly detected, suggesting that the collapsing core is formed locally inside a 
flattened structure (§4.2). 

4. The exhibited hourglass-like B^ morphologies i n the e2 and possib l y the e8 dust ridges 



are consistent with the propos ed local collapse (|Ho &: Young 19961 ; ISollins et al. 2004 ) 



instead of the global collapse (IRudolph et al. 19901 ) . This indicates that both the e2 



and e8 cores are in a supercritical stage (gravity dominating over the B field) at the 
0.02 pc scale. In contrast, the B field morphology of the 0.5 pc envelope is uniform 
across the e2 and e8 dust ridges and strong (>1 mG). We further propose that the 
B field in the 0.5 pc scale envelope is subcritical (B field dominating the gravity). 
That the supercritical cores forme d inside a subcritical envelope seems to support 



a magnetic frag mentation scenario (IMouschovias &: Morton 19911 ; iLizano &; Shu 1989 



P 



Shu et al. 2004 ). suggesting that ambipolar diffusion plays a key role in the evolution 



of the envelope at this stage (§4.5). 
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Appendix 

In order to test how significantly the measured B field deviates from a purely radial field, 
we have analyzed the differences, 5, between position angles (P.A.s) of the measured B field 
lines and the corresponding hypothetical radial field lines. The radial field lines are derived 
from their relative positions (center of each P. A.) with respect to the origin, which is defined 
as the 870 /iin peak in e2. Due to the limited detected data points, we do not apply this 
statistical analysis to e8. All the data considered are above 3 07 . We have excluded the 6 
data points in the north-west extension, which correspond to e2 pol NW. The distribution of 5 
apparently deviates from a Gaussian (Figure Al (Left-panel)). When, nevertheless enforcing 
a Gaussian fitting, the derived mean ji is 1.6° and the standard deviation a is 14°, which is 
larger than the mean measurement uncertainty <7 me an,PA of 6.9°. 

To quantify the deviation from a Gaussian, we further apply a Kolmogorov-Smirnov 
(KS) test to 5. As null hypothesis we assume that the distribution of S is normal with <7 mea n,PA 
and /i — 0°. This mimics an observation of a radial field with our observed measurement 
uncertainty o" mca n,PA- The measurement uncertainties in each PA. are propagated with 
a Monte-Carlo (MC) simulation when deriving 5 and applying the KS test. In the MC 
simulation we also allow for a shift of iO'/l of the origin of the radial field. Figure Al 
(Right-panel) shows the cumulative distributions which are used for the statistic measure in 
the KS test. As a result, the probability of the measured field being radial is 20%. For the 
ilc"mcan,PA error bounds, the probability of being radial is less than 5%. Therefore, the null 
hypothesis can be rejected. Our test favors the existence of a non-radial B field. 

5 as a function of PA. is presented in Figure A2. If we separate the segments according 
to the depolarization zones (marked as arrows in Figure A2), there are always segments with 
S larger than 3 o" mcan pA in each zone. This suggests that deviations are not prevalent in 
certain directions, but rather grouped together and interleaved with depolarization zones, as 
we would expect from an hourglass field morphology. 
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Table I. SMA dust polarization at 870 jum in e2 
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Table 1 — Continued 



Ax 


Ay 


I 


P{%) 




0(°) 


0.3 


0.0 


0.60 


2.5 ± 0.6 


15 


80 ± 7 


0.0 


0.0 


0.95 


2.1 ± 0.4 


20 


-79 ± 5 


-0.3 


0.0 


0.76 


1.7 ± 0.5 


13 


-80 ± 9 


0.0 


-0.3 


0.23 


6.0 ± 1.6 


14 


-81 ± 8 


-3.6 


-2.4 


0.22 


5.0 ± 1.7 


11 


-17 ± 10 



Note. — Ax & Ay: offsets in arcsecond from Right 
Ascension (J2000) = 19 h 23 m 43 s .95, Declination (J2000) 
= 14°30'34'/00. I: intensity of the Stokes I component 
in Jy beam . P(%): polarization percentage, derived 
from the ratio I p /L I p : the polarized intensity in mJy 
beam -1 . <fi: position angle from the north to the east. 
Listed data points are within the red box associated 
with e2 in Figure 2. All listed points are above 3 ctj . 
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Table 2. SMA dust polarization at 870/xm in e8 
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2.2 ± 0.6 


14 


-1 ± 8 


-1.5 


-5.7 


0.22 


5.8 ± 1.8 


13 


-1 ± 9 


-0.9 


-6.0 


1.57 


0.7 ± 0.2 


11 


-70 ± 10 


-1.5 


-6.6 


0.27 


7.5 ± 1.4 


20 


-21 ± 5 


-1.5 


-6.9 


0.39 


3.8 ± 1.0 


15 


-19 ± 7 


-1.2 


-7.2 


0.55 


2.0 ± 0.7 


11 


36 ± 10 


-1.2 


-7.5 


0.42 


3.3 ± 0.9 


14 


59 ± 8 


1.5 


-8.7 


0.26 


4.2 ± 1.4 


11 


6 ± 10 


-1.8 


-8.7 


0.24 


5.0 ± 1.6 


12 


-33 ± 9 


-2.1 


-8.7 


0.23 


8.0 ± 1.7 


18 


-40 ± 6 


-2.1 


-9.0 


0.23 


4.7 ± 1.5 


11 


-33 ± 10 



Note. — Identical notation as in Table 1. Ax 
and Ay are with respect to the same coordinate 
as in Table 1. Listed data points are all within 
the red box associated with e8 in Figure 2. 



Table 3. Derived parameters in e2 and e8 



Parameters units e2 e8 





BIMA a 


SMA 6 


BIMA a 


SMA fe 




n 


3.6 x 2.6 


0.9 x 0.8 


5.0 x 2.2 


0.9 x 0.3 


Mr 


M 




220 




94 


n H2 


cm -3 


1.5 x 10 6 


2.7 x 10 7 


7.6 x 10 5 


1.2 x 10 7 




mG 


>1 


<19 


>1 


<8 


Aa 


mpc 


>5 


<6 


>H 


<5 




(") 


(>0.2) 


(<0.2) 


(>0.3) 


(<0.2) 


At,ct 


mpc 


25 


6 


36 


9 




(") 


(0.8) 


(0.2) 


(1.1) 


(0.3) 


Am,ct 


mpc 


>24 


<25 


>47 


<24 




(") 


(>0.7) 


(<0.8) 


(>1.4) 


(<0.7) 



"The derived mean nn 2 is based on a Gaussian fit of the 1.3 mm continuum emission toward the peak. All 

the natural length scales depend on this nn 2 • 
b The derived mean nn 2 is within a radius of 1" centered on the 870 /xm continuum peak, where the 
collapsing signatures were clearly revealed with the molecular lines. All the natural length scales and the 

upper limit of B depend on this «h 2 ■ 



Note. — dmaj and Q m %n refer to the major and minor axis of the deconvolved 
size, respectively. They are determined by a best-fit Gaussian of the continuum 
emission centered on the peak. M R refers to the estimated M gas from the contin- 
uum emission within a radius of 1" centered on the peak position. The lower limit 
of the B field strength of 1 mG is adopted from Lai et al. (2001). Characteristic 
length scales Aa, At,ct and Am, ct are calculated assuming T = 100 K. Note that 
the weak constraints on Aa and Am,ct result from the large possible range of B field 
strength. 
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Fig. 1. — Schematic cartoon of structures in W51 e2/e8. e2 and e8 refer to the UCHII 
region e2 and e8, respectively. 
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(a) BIMA 1.3 mm W51 e2/e8 (b) SMA 870/rni W51 e2/e8 




offset (arcsecond) offset (arcsecond) 



Fig. 2. — (a) Polarization map observed with BIMA adopted from Lai ct al. (2001). Contours are the 
1.3 mm continuum strength of —3, 3, 5, 10, 20, 30, 90, 100, 110 x 27 mJy beam -1 , where the size of 
the synthesized beam is 2'.'7 x 2'.'0. The cyan stars mark the UCHII regions. The offset is in arcsecond with 
respect to Right Ascension (J2000)=19 h '23 m 43 s .95, Declination (J2000)=14°30'34"00. Red boxes indicate 
the regions presented in panel (c) and (d). (b): SMA 870 fim continuum emission in W51 e2/e8. Contours 
are the 870 /xm continuum strength at —6, —3, 3, 6, 12, 24, 36, 48, 60 x 60 mJy beam -1 , where the beam 
size is Of! 7 x 0'.'6. All the other symbols are identical to the ones in (a), (c): SMA polarization in W51 
e2. Contours are the 870 /im continuum strength at —4, —2, 2, 4, 6, 8, 10, 15, 20, 55, 60 x 60 mJy 
beam -1 . Black and red segments represent the polarization with its length proportional to the polarized 
perce ntage at 2 to 3 aj r and above 3oj p , respective ly. Solid and unfille d triangles mark the positions of the 
H 2 (jGenzel et al. 198 lb and (J,K) = (9,6) NH 3 |Pratap et al. 199lh masers, respectively. All the other 
symbols are identical to the ones in (a), (d) SMA polarization map in W51 e8. All the symbols and contour 
levels are identical to the ones in (c). All images shown have been corrected for the effect of primary beam 
attenuation. In each panel, the synthesized beam is plotted as a black ellipse at the lower-left corner. In 
panel (a), (c) and (d), the polarized intensity is shown in color scale with strength indicated by the wedge 
on the right in units of mJy beam -1 (mJy/B). 
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8 0-8 
offset (arcsecond) 

Fig. 3. — Polarization map of the SMA results (red segments) at 870 /im restored with the 
BIMA synthesized beam. BIMA polarization at 1.3 mm is shown in blue segments. The 
synthesized beam size is 2'. / 7x2'.'0 with a P. A. of 1°. The black contours represent the SMA 
870 /im continuum emission after restoring to the BIMA resolution with contours of —3, 3, 5, 
10, 15, 20, 25 x 0.3 Jy beam -1 . The grey contours represent the BIMA 1.3 mm continuum 
emission with strengths as shown in Figure 2. All polarization vectors plotted are above 
3cr/ p . All the other notations are identical to the ones in Figure 2. The slight offset between 
the SMA and BIMA maps is most likely due to missing short spacing information in the 
SMA data, being insensitive to structures larger than 8". 
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Fig. 4. — Plot of polarization percentage, P(%), versus normalized intensity (I /I max )- The 
blue circles, red crosses, and black triangles are SMA results for e2, e8 and BIMA results for 
both e2 and e8, respectively. All the plotted data are above 3a i. 
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19 h 23 m 44 s 43 s 9 43 s 8 19 h 23 m 44 s 43 s 8 

a(J2000) a(J2000) 



Fig. 5. — Left-panel: B field maps in W51 e2 with S/N ratio between 2 to 3 aj p (black 
segments) and >3 o/ p (red segments), and 870 /im continuum emission in black contours. The 
dot-dash cyan segment marks the direction of the ionized accreting flow by Keto & Klaassen 
(2008). Blue segments are the inferred B field direction from the BIMA. The magnetic field 
segments (all of equal length) are rotated by 90° with respect to the polarization segments. 
Right-panel: B field map in W51 e8. All the symbols are identical to the ones in the 
left-panel. All the contours are plotted with the same level as in Figure 2(c) and 2(d) for 
left-panel and right-panel, respectively. Note that even at 2 to 3 cr/ p , the B field directions 
are varying in a coherent manner. Hence the trends are significant. 
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t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




19 h 23 m 44 s 5 44fO 43f5 

a(J2000) 

Fig. 6. — B field maps of W51 e2/e8 from the SMA (red segments) at 870 /im and BIMA 
(blue segments) at 1.3 mm. The black and grey contours represent the strengths of the 
continuum emission at 870 /im and at 1.3 mm, respectively. The black contours are plotted 
in the same levels as in Figure 2(b), and grey contours plotted are 3, 5, 10 x 27 mJy beam -1 . 
The other symbols are identical to the ones in Figure 2. The synthesized beams of the SMA 
and BIMA are plotted in the lower-left corner as white and black ellipses, respectively. 
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E 




S (°) S (°) 

Fig. Al. — Left-panel: Histogram of the differences 5 between the measured P.A. and a 
radial field. The cyan dashed curve is the enforced fitted Gaussian to the measurement with 
the resultant standard deviation a=14° and mean /i=1.6°. The red solid curve is the normal 
distribution with a = cr meani pA=6.9°, which is the mean uncertainty of the measurement. 
Right-panel: Cumulative distribution of 5 (solid blue line). The red curve marks the cu- 
mulative normal distribution with 0" meani pA- The ± 1 cr mca n,PA error bounds from the MC 
simulation, including measurement uncertainties and the shift of the center of the radial 
field, are also shown. 
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Fig. A2. — Plot of position angle versus 5. The red lines mark the 3 0" meani pA bounds. 
The black arrows mark the depolarization zones around P.A.s of 60° and 220°, ~ parallel 
to the maximum velocity gradient in the H53a line. The areas with respect to the 870 
continuum peak are marked as NW, NE, SE and SW, which correspond to the north-west, 
north-east, south-east and south-west, respectively. There are always 8 > 3 o~ moa n,PA m the 
areas separated by the depolarization zones. 



